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Abstract. - The odd-parity spin-triplet Cooper pairing states show a six-fold degeneracy for 
a quasi- two-dimensional electron system. In this article we show how spin-orbit coupling can 
lift this degeneracy for Sr2Ru04 taking the band structure based on the three relevant t2g- 
orbitals of the Ru-ions into account. The stabilized state depends on the relative strength of 
the pairing interaction. We show that under reasonable assumptions the chiral pairing state 
d(k) = z{kx ± iky) is favored against the others. 



The search for new superconducting transition metal oxide systems has lead a few years 
ago to the discovery of Sr2Ru04 Despite the low transition temperature, {Tcmax ~ 1.5K ) 
which is, in addition, very sensitive to disorder |^ , the unusual properties of the superconductor 
have attracted much interest. During the last year it has been experimentally established 
that this compound is a spin-triplet superconductor as initially suggested on theoretical basis 
[ 3 , Q . The superconducting state is characterized by the violation of time-reversal symmetry 
[5| and equal spin pairing within the basal plane of the tetragonal crystal lattice The 
single candidate consistent with all presently known experimental data is described by the gap 
function of the symmetry d(k) = z(kx ± iky) in the standard vector notation This state is 
two-fold degenerate and has chiral symmetry, i.e. the Cooper pairs possess an orbital angular 
momentum. There are various possible other candidates for spin-triplet pairing which might 
be rather close in energy as we will show below. The question arises what is the microscopic 
reason to favor the chiral phase compared to the others. It is the aim of this paper to show 
that the chiral state is a natural consequence of the electronic structure of Sr2Ru04 and the 
presence of spin-orbit coupling. 

Sr2Ru04 is like several of the high-temperature superconductors a layered perovskite system 
of stacking Ru02-planes with a single layer per unit cell [Q. Like Cu in the high-temperature 
superconductors Ru forms a square lattice. While a single-band description seems to be 
adequate for the Cu02-system the case of Ru02 requires at least three bands originating from 
the three 4d-i2g-orbitals of Ru occupied by four electrons on the average. All bands cross the 
Fermi energy and give rise to three Fermi surfaces, two electron-like and one hole-like, all of 
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nearly cylindrical shape owing to the weak dispersion along the c-axis ]9|, |10|, 11 . The low- 
temperature normal state of Sr2Ru04 is that of a strongly correlated, nearly two-dimensional 
Fermi liquid. In this sense this system resembles to some extent the nature of the Fermi liquid 
'Hei. 

Starting from this assumption we show that the ordinary weak coupling theory for spin- 
triplet pairing yields a six-fold degenerate superconducting state if the spin space is completely 
rotation symmetric. This degeneracy is lifted through spin-orbit coupling as suggested in 
Ref.(|^, H). In this letter we describe the effect of spin-orbit coupling based on the orbital 
origin of the three electron bands. For this purpose we introduce a simple model taking the 
basic structure of the three orbitals into account and leading to spin-triplet superconductivity. 
The three i23-orbitals, dyz, dzx, dxy, will be labeled as 1,2,3, respectively. The kinetic energy 
is then given by 



Hfe™ = -^Ct I ffk ek2 Ck. (1) 

k,. \ £k3 / 

with = {cliscli2s^k3s)^ (eki,ek2) = - 2ti {cos ky , cos k^) - £k3 = -2t2icoskx + cosky) - 
Ats cos kx cos ky — jj! and ^k = —4^4 sin kx sin ky (TV: number of lattices sites]. The parameters 
will be later chosen to fit the band structure as close as possible ||l^, |l^. Through the 
hybridization of the orbital 1 and 2 the Hamiltonian leads to the bands a and (3 corresponding 
to the hole-like and electron-like Fermi surface, respectively. The band belonging to orbital 3 
is called 7 in accordance with the band structure literature |lO[ |l3j . 

We now turn to the discussion of the pairing interaction yielding spin-triplet pairing. The 
onsite interaction is definitely repulsive suppressing pairing in the standard s-wave channel. 
Therefore we focus on nearest-neighbor-site interactions. Furthermore, the zero-momentum 
pairing requires that two paired electrons are on the same Fermi surface. Interband pairs have, 
in general, a finite net momentum. As a possible candidate for pairing ferromagnetic spin 
fluctuations were discussed [|[ @]. Experimental studies show, however, that the ferromagnetic 
component of the spin fluctuations are rather weak [T5| and, in additions, if they are enhanced, 
e.g. by pressure, superconductivity is actually suppressed ||l^. Alternatively, Hund's rule 
coupling was mentioned as another cause for parallel spin pairing ||, ^. Standard Hund's rule 
coupling on the Ru-ion, however, yields an onsite interaction and would not be available for 
spin-triplet pairing. There is, however, a spin-dependent interaction mediated via orthogonal 
oxygen orbitals between different nearest-neighbor orbitals. For example, on a bond between 
site 1 and 2 along the x-axis spins on the Yln-dxy- and da:2-orbital on different sites interact 
ferromagnetically as they hybridize with orthogonal p-orbitals on the intermediate oxygens. 
This is a variant of the Goodenough-Kanamori mechanism for ferromagnetic exchange based 
on Hund's rule effect on oxygen in transition metal oxides. It is easy to see that the d^^y-orbital 
couples in this way to the dzx- and dyz-orbital of the Ru-ions on neighboring sites along the 
X- and y-direction. 



Hint — J' / ^[Si3 • Si+x_2 + Si+x,3 ■ Si2 + Si3 • S,;+y^i -|- S^+y^a • S^i] (2) 

i 

This is an interorbital interaction and by itself for pairing not so effective. However, the local 
spin polarization generated by this interaction can be transfered via onsite Hund's rule coupling 
to the other orbitals and lead to an effective spin dependent nearest neighbor intraorbital 
pairing interaction which has the form 
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Hp — ^^[Ji{Sii ■ Si+xi +Si2 •Si-|_y2} + '^2{Sil • Si+yi +Si2 •Si-|_x2} + >/3{Si3 -Si+xS + Sjs • Si+ys}] 
i 

_ _ (3) 

where Ji is ignored as it has no contribution in our modeL The effective couphng strengths 
depend on the onsite polarizabihty of the spins in the different orbitals. 

./2 = -^^^51^33 (q) and J3 = -^^-^^[Xii{fi) +Xi2{tl)] (4) 
q q 

where Xi-f^il) = Jq^ dt{[Sf^ci{t),S^^q{0)]) is the static spin susceptibihty (ignoring retar- 
dation effects) of the orbital /i with the field acting on the orbital ly. 

Besides this spin dependent contribution we introduce also a spin independent part of the 
attractive intraorbital interaction of unknown origin. This can be done by replacing S^^ • Sj^ 
by Si^ ■ Sjp + v^Uiij^nj^/A. Decomposed into spin-singlet and spin-triplet channel, this does 
not change the structure of the interaction. For simplicity, we will take vi = V2 = W3 for 
our calculation. To lift the degeneracy completely we have also to include the pair scattering 
between different Fermi surfaces. This interaction is mediated by Coulomb interaction [ p^ . 
For unconventional Cooper pairs the onsite scattering is ineffective. Thus we include the 
scattering via nearest-neighbor interaction also here, 

for n ^ fj,' . Including all terms and changing to momentum space we obtain in a pairing 
interaction of the form, 

Hpair — ^^i',si,...S4 ~ ^r)Ck^siCLk^S2''~k>'S3CkVs4 (6) 

where indicates that we sum over all repeated indices (r — x,y) and the coefficients 
^/T/i' si S4 are easily obtained from the above interactions. The decomposition of cos(fc — k') 
into cosfccosfc' -I- sinfcsinfc' yields the wave functions of the possible states: "extended s- 
wave" (cos kx ~t- cosfcy), "d2;2_y2-wave ' (cos/c^ cos/cy) for the spin-singlet channel and the 
p-wave" spin-triplet states (sin fc^,, sin fcj^) which are degenerate by symmetry While both 
spin-singlet states have, in general, zero-nodes in the gap, it is possible to combine the two 
degenerate p-wave components to form a nodeless gap. In this sense the p-wave channel is 
expected to be favored, if the spin dependent part of the interaction does not suppress the 
triplet spin configuration. 

The degeneracy of the two Cooper pair channels with the orbital wave functions sin and 
sin ky is independent of the spin wave function. Thus, we may choose the spin configuration 
to optimize the condensation energy. The gap function is given by 

A^(k) = io- • d^(k)(T^ = i ^ CT'^a^i>(a,,^j sin fc^^ -I- 6^^ sin fcj^) (7) 

u—XjyjZ 

where P is the unit vector in i^-direction. The quasiparticle energy gap in each band is given 
by v^|d^(k)|2. States with the same gap have within the weak coupling theory the same 
condensation energy and are consequently degenerate. The combinations given in Table I 
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Table I. - Degenerate stable spin-triplet pairing states listed according to the irreducible representa- 
tions of the tetragonal point group D^h assuming simultaneous rotations of spin and orbital part. 



r 




P Q 


Am 


X sin fca; + y sin ky 


+1 +1 


Am 


X sin ky — y sin k^ 


+1 -1 


Biu 


X sin kx — y sin ky 


-1 +1 




X sin ky y sin kx 


-1 -1 


Eu 


z(sinA;2; ± i sin ky) 


±1 



have all the same gap y sin^ kx + sin^ ky so that we find six degenerate pairing states ^ . 
This is a special feature of a two-dimensional Fermi liquid. In three dimensions one finds only 
one stable state, the Balian-Werthamer state, d(k) — xs'mkx + ysinfcy + zsinfc^. Note, 
that this result is general and does not depend on the details of the Hamiltonian as long as 
it favors p-wave pairing. It was shown that going beyond weak coupling by including the 
renormalization of the spin fluctuation spectrum the state d(k) = z{smkx ± i sin ky) would 
be favored analogous to the spin fluctuation feedback mechanism stabilizing the A-phase of 
superfluid '^He This is a secondary effect and does not affect the degeneracy of Tc, but is 
connected with the condensation energy at low temperature. 

The degeneracy of Tc can be lifted by removal of spin rotation symmetry, i.e. by taking 
spin-orbit coupling into account. The spin-orbit coupling is connected with the heaviest ion - 
Ru in our case - where it acts as an onsite term in the Hamiltonian, Tiso ~ -^X^i ' ^i- The 
angular momentum operates on the three t2g-orbitals on the site i. If restrict ourselves to 
these three orbitals, ignoring the e^-orbitals, we can write the spin-orbit Hamiltonian as 



l.rn,n W,s^s' 

where eimn is the completely antisymmetric tensor and A is a coupling constant such that the 
t2£,-states behave like an ^ = 1 angular momentum representation (we use for /, m, n either 
{x,y,z} or equivalently {1,2,3}). Together with Eq.(^ the spin-orbit coupling leads to new 
quasiparticles which are now labeled by pseudo-orbital and pseudo-spin indices connected with 
the original ones by the unitary transformation defined for each wave vector k, 



(Okls, ak2s, flkSs)^ = (Ckls, Ck2s, Ck3-s)^C^ks (9) 

where the spin-orbit coupling mixes the up spins of the orbital 1 and 2 with down spin of 
orbital 3 and vice versa. We find that the matrix elements of C/ks of opposite spin are related 
as, Uk/x^s = ul_^^ _^ for /i = 1,2 and Wksi^s = -""ksi^.-s ^ = 3. Note, that the 

inclusion of spin-orbit coupling affects the band structure so that the chemical potential of 
the orbitals have to be adjusted for different values of A to produce the proper Fermi surface 
shape. 

Cooper pairing has now to be defined and analyzed concerning the symmetry based on these 
new quasiparticles. We proceed, therefore, by rewriting the original pairing interaction Hpair 
in terms of these quasiparticles, 
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Hpair — ^si's2S3S4('*^''* )^ki/si'^-ki^S2'^-''''^'*3'^kVs4 (10) 

where 

siS2S3S4*^ ' / I -'^81828384 '^kmi^siS2 k'm'i/'s3S4 ' -'^,siS2S384'^km;/si82 k'm'i/'s384 y 

-^81828384 ~ Jm5mm"^siSi'^S2S3 ~^ 9mm' ^8184^8283 + ffmm' ^81 ,— 84 <^8i ,— 82 '^83, — 84 

ij-mm' = a' ,S 8 8 Cll'l 

81828384 i/mm' "81, -84 "81 ,82 "83, 84 K'-'-l 

Here we introduced the abbreviation J„i = (J„i/4)(1 + v). g and g' are symmetric with 
non-vanishing elements gi2 = g'l^ = g'23 = g- The functions , P = ±1, are used to express 
the pairing states and shall be constructed according to the symmetries listed in Table I where 
P for each state is given, 

/P 



^kli.8182 = ■"klysiMkliys2 Sinfcx + «^Mk2i.siWk2i.S2 sin/Cy 

^k2i/siS2 = Mk2i.siUk27/s2 sin fcx + iPwkiiysi Mkiiys2 sm ky 

^k3i^siS2 ^ Mk3i/siMk3,y82(sinfcj; - isiS2Psinfcj,) (12) 



A 90° rotation acting on the orbital part only (i.e. k) leads to l^kmvsiS2 ~^ ~*-^^kmi/si82 
following Table I. Note, that in Hpair different signs of P are coupled through the pair scattering 
between orbital 3 and orbital 1 or 2 only. These intcrband pair scatterings are responsible for 
the lifting of degeneracy of the Ai^2u and -Bi,2m pairing states. 

The standard BCS-type mean-field approach leads to the definition of the quasiparticle gap, 

A.828l(k) = Vr:l,,.3,.,(k,k')/.'8384(k') (13) 

with /j^/s3s^(k') = (a_k'i/'83ak'j/'84) which is readily calculated after diagonalizing the mean- 
field Hamiltonian. After some algebra, we formulate the self-consistent gap equation for pairing 
in the odd-parity (pseudo spin-triplet) channel for the d-vector. 



= {jmSmm' + ~9mm' ' Q~9'mm')) {QU^mtn + "k^.Tt) 

'^^(k) = ^]^E {jm^mm' + 9mm' — Q9mm'^ (^~Q^kmull +^krnj/TT) 
dtO^) = (j^m^mm' + 9mm' + g'mm'^ ^kmul l^m' (14) 

with and D^, are, respectively, the eigenvectors of 

-^m? = ~]yE i^^m^mm" +gmm" — Qg'mm')^^m'v'-[-[^^mv''\-[Pv'i!.D^, 
^m' = ~]yE ij^™-^™-"'-'' ~^ 9mm'' + g'mm'^^^7n'v'^i^^mv'^iPi''^D'i^n (15) 

where Fv,\s. = tanh(/3i?i,/k/2)/2£'y/k and P,Q = ±1 corresponding to Table I. These six 
eigen-equations (4 for and 2 for D^,) has solutions with the symmetries corresponding 
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to the pairing states as in Table I. Note that Df^, is independent of Q and is the same for both 
P = ±1 (because Wj^^j^-f^ — l^urniyU '^^i^^ therefore, always doubly degenerate. If D^, 
gives the highest transition temperature Tc after solving the eigen-equations self-consistently, 
then all has only trivial solutions at Tc- Hence the pairing state has only component 
and is doubly degenerate, which is consistent to the chiral state. On the other hand, there will 
be no component in the pairing state when one of the gives the highest Tc. 




Fig. 1. - Phase diagram of stable states, with phase I represents z{kx ±iky), phase II Hky — ykx and 
phase III yik^+yky. The parameters used are U = 0.18 eV, t2 = 0.17 eV, ts = 0.08 eV, ti = 0.02 eV, 
fi' - fj. = 0.1 eV, Js = -0.1 eV and g = 0.02 eV. 



Our calculation shows that the spin-orbit coupling lifts the degeneracy in favor of different 
pairing state depending on the relative strength of the pairing interactions, J2 and J3. In Fig.l 
we present the phase diagram of spin-orbit coupling A versus the ratio J2/ Js- Among the three 
states occurring in the phase diagram we find also the time reversal symmetry breaking state 
z(kx ± iky) which is most stable, if the pairing interaction is dominant in the d^y-orbital. 
The other two phase corresponding the irreducible representations yli^ (phase II) and 
(phase I) in Table I, require either substantially stronger pairing interaction in the bands of 
the orbitals dyz and d^x or rather weak spin-orbit coupling. The values of A of the Ru^+-ion 
in the literature are of the order of 0.07 eV. In view of the effect the spin-orbit coupling 
has on the quasiparticle spectrum we may argue that its strength to be moderate, since the 
band structure calculations ignoring spin-orbit coupling are in reasonable agreement with the 
de Haas- van Alphen data |l^, There are also experimental indications which suggest 
that superconductivity in Sr2Ru04 is carried mainly by the dxy-hand Jl^ . Clearly also the 
electronic density of states of this band makes up 43% of the whole system. (In the choice of 
parameters used in this calculation we have taken the ratios of density of states into account.) 
The role of the interband pair scattering in Eq.(^ is to distinguish between the Ai^- and 
^2«-state which would be degenerate for 5 = 0. 

Our approach is based on the assumption that spin-orbit coupling is not too strong so that 
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we can deal with it essentially by changing only the character of the quasiparticlcs. This 
procedure is complicated due to the itinerant character of the electrons. In our discussion 
we have not taken into account the modification of the spin dependent matrix elements of 
the interaction. For this purpose we would need more detailed knowledge about the pairing 
mechanism. Nevertheless, for the mechanism based on the extended Hund's rule coupling we 
can test the effect of spin orbit coupling. Calculating the free-quasiparticle spin susceptibilities 
including spin-orbit coupling shows that the anisotropy in the spin interaction in Eq. is only 
of the order of 1% in the range of A important for, considerably smaller than the above effect. 

Calculating the static spin susceptibility for all bands we obtain within our model a very 
similar form as calculated from LDA and measured recently by neutron scattering [p^ . 
Both theory and experiment agree well in the aspect that the strongest spin correlation is not 
ferromagnetic, g = 0, but has a finite wave vector associated with the rather strong nesting 
feature of the nearly one-dimensional bands of the dyz- and dzjj-orbitals which in this way 
dominate the g-dependence of the susceptibility. The large density of states, however, leads 
to the strongest contribution to the overall susceptibility by the dxy-h&nd. In calculating 
coupling strengths J2 and J3 we find that their ratio is of the order 1.7 in favor of J2 which 
still would lead to the chiral superconducting state in the phase diagram. However, also the 
spin independent interactions have to be taken into account, which we cannot estimate at 
present. The main conclusion of our study is that spin-orbit coupling plays an important role 
in stabilizing the time reversal symmetry breaking p-wave pairing state in Sr2Ru04 and is 
strongly connected with the orbital structure of the electron bands. 
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